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The present work analyzes the importance of the different components of the environment effects 
on the local spectroscopy of extended strongly correlated systems. It has been found that the 
usual formal charge definition of the charge transfer and Madelung potential are far too crude for 
an accurate determination of the local excitation energies in embedded fragment calculations. A 
criterion for the validation of the embedding against the infinite system density of states has been 
proposed. 



INTRODUCTION 



In the last twenty years the science of strongly cor- 
related materials has attracted a lot of attention. In- 
deed, chemists have synthesized a large number of new 
families of materials which present unusual and fascinat- 
ing properties directly related to the strongly correlated 
character of their electronic structure (to cite only a few 
systems, one can think of the copper oxides presenting 
d-wave, high-Tc super-conductivity 0, manganites with 
giant magneto-resistance [31, Prussian Blue analogs with 
photo-magnetic properties [2j, etc.). In all these materi- 
als the electrons responsible for their spectacular prop- 
erties are few (per unit cell), usually unpaired and lo- 
calized both spatially and energetically near the Fermi 
level. That is to say that the wave function of these sys- 
tems is essentially multi-configurational and cannot be 
correctly treated by single-reference based methods such 
as Hartree-Fock plus perturbation theory or even Density 
Functional the only ah initio methods to be tractable 
on infinite crystals. Physicists have thus relied on semi- 
empirical models such as the Heisenberg t — J 
or Hubbard Hamiltonians in order to describe these 
strongly correlated materials. These effective models are 
of valence-bond (VB) type and belong to the zero differ- 
ential overlap (ZDO) class. They are aimed at describing 
the effective microscopic interactions between the set of 
valence electrons responsible for the macroscopic low en- 
ergy properties. For instance, this is in this respect that 
the famous t — J model has first been derived by Zhang 
and Rice et al ^ a,s a model for the high-Tc copper ox- 
ides. Indeed, the t — J model describes the dominant 
interactions between the magnetic electrons of the Gu02 
superconducting planes. These planes can be seen as a 
Qi^+ square lattice, each copper ion supporting (before 
doping) one dx2_y2 magnetic electron. These unpaired, 
Fermi level electrons (considered to be the important 
ones for the superconductivity properties) only interact 
between each others through super-exchange processes 
(via the bridging oxygen anions 0^~) and through effec- 
tive hopping integrals (when the planes are doped as in 



the super-conducing phases), i.e. through a i — J type of 
Hamiltonian. 

One sees immediately that such representations neces- 
sitate i) to clearly identify the pertinent electrons and 
orbitals responsible for the macroscopic low energy prop- 
erties — they will be called from now on active electrons 
and active orbitals — ii) to determine the dominant in- 
teractions between these electrons iii) finally to evaluate 
their relative amplitudes. While the nature of the active 
electrons and orbitals, and even the form of the domi- 
nant interactions, are most of the time relatively easy to 
derived using simple chemical considerations, the quan- 
titative determination of the interactions amplitudes is 
a much more troublesome problem, even-though as cru- 
cial. Indeed, the low energy properties of the macroscopic 
systems are highly dependent on the relative amplitudes 
of the different effective integrals, since the latter can 
govern properties as important as the metallic versus in- 
sulating character of a compound. Unfortunately, ex- 
perimental data are, most of the time, unable to totally 
determine them and it is necessary to rely on ah initio 
quantum-chemical calculations for this purpose. Since 
these effective interactions (hopping, magnetic exchange, 
on-site bi-electronic repulsion, etc.) are essentially local 
(see references 8, 9] for detailed analyses of the locality 
of the different effective interactions) and involve a small 
number of bodies (most of the time only one or two), 
they can be accurately determined from the computa- 
tion of the local excitation energies and related states. 
For instance, in the high-Tc copper oxides, the magnetic 
coupling between two active dx^^-yi electrons can be de- 
termined as the local singlet to triplet excitation energy 
between two nearest neighbors, oxygen bridged, copper 
atoms. The embedded cluster technique, where the local 
valence excitations are computed on an adequate frag- 
ment of the crystal, has proved to be one of the most ef- 
ficient technique for this purpose. Indeed, the full power 
of the ah initio molecular-spectroscopy methods can be 
used on such finite systems. One can cite the successes 
obtained on high-Tc copper oxides vanadates ox- 

ides 0, 0|, fluoroperovskites 0, etc. For instance, 
the magnetic exchange in the high-Tc superconductor 
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Nd2Cu04 has been computed within experimental accu- 
racy (Jca/c, = — 126.4 mey being the computed value in 
reference [iflj and Jexp. = —126 ± 5meV being the mea- 
sured experimental value [l^'). 

As we will see in details later on, the local character of 
the effective interactions refers to the local character of 
the electronic processes involved in the local excitations, 
once the mono-electronic part of the wave function is de- 
fined, that is, once the shape of the fragment orbitals 
and their energies are defined. It is however obvious that 
the latter, and in consequence the effective interactions 
amplitudes, are strongly influenced by the fragment crys- 
talline surrounding. Indeed, the orbitals and orbital en- 
ergies of an isolated fragment are very different from the 
orbitals and orbital energies of the same fragment in a 
crystalline environment. It is therefore crucial to embed 
the studied fragment into an environment reproducing 
the main effects of the rest of the crystal. 

The purpose of this paper is to analyze, review and 
clarify the requirements for a good embedding and in 
this light to discuss the most common practices. We will 
also propose a general technique for checking the quality 
of an embedding and study the sensitivity of some ef- 
fective parameters, such as the magnetic super-exchange 
and the effective hopping integrals, to the quality of the 
embedding and the crystalline environment. The next 
section will therefore be devoted to the analysis of the 
environment-fragment interactions and the proposal of a 
simple criterion measuring the embedding quality. Sec- 
tion three will review the usual practices, exemplify the 
influence of the crystalline environment on real systems 
and point out the larger importance of the different em- 
bedding components on the effective parameters than is 
commonly assumed. Finally, the last section will be de- 
voted to conclusions and perspectives. 

INTERACTIONS BETWEEN THE FRAGMENT 
AND THE REST OF THE CRYSTAL 

The interactions between a selected fragment and the 
rest of the crystal can be classified into two different 
types : the effects on the fragment of the crystal average 
electrostatic field and the effects of its fluctuations. 

Let us first start with the second kind : the electro- 
static and spin fields fluctuations, that is the instanta- 
neous response of the rest of the crystal to a particular 
electronic configuration of the fragment. In other words 
we are talking about the dynamical correlation effects of 
the rest of the crystal on the fragment valence states un- 
der consideration. The spatial separation between the 
active centers and their crystalline environment allows 
us to conduct a multi-polar analysis of the different con- 
tributions to the excitation energies, following the same 
pattern used in the description of the inter-molecular in- 
teractions. Along this line, it is very simple to see that 



the leading dynamical-correlation contributions are the 
polarization responses of the environment first to a local 
charge, then to a dipole moment or transition dipole mo- 
ment on the fragment In other words, these effects 
correspond to excitations from the average environment 
(chosen as its zeroth order description) in response to 
a particular fragment configuration. In a perturbative 
scheme they can be described as the coupling between 
single excitations in the environment to a local charge, 
a local dipole or a local transition dipole on the frag- 
ment. From multi-polar expansion combined to second 
order perturbation theory, it is easy to see that the lead- 
ing contributions come from the interactions between a 
charge on the fragment and a dipole on the environment. 
Such effects decrease as (1/i?^)^, that is quite rapidly. 
For a complete description and analysis of the differ- 
ent contributions one can refer to Provided that 
the largest l/R^ terms — i.e. the terms coming from 
single-excitations on the first neighbors of the magnetic 
centers — are correctly treated, the rapid decrease of 
these effects as a function of R allows us to consider the 
rest of them as reasonably negligible. Indeed, while the 
first neighbor shell of atoms surrounding the magnetic 
sites are usually at a distance of 2A to 3A(from 3.8a.u. 
to 5.7a.u.) the second neighbor distances are usually of 
the order of SA to SA, (5.7a.u. to 9.5a.u.). The sec- 
ond neighbor contributions can therefore be expected to 
be of an order of magnitude 10^^ to 10~^ smaller than 
the on-site polarization contributions. One can therefore 
consider that including in the studied fragment, the first 
coordination shell of all the atoms involved in the local 
excitations under consideration, is enough to take into ac- 
count the essential contributions to the modification of 
the effective interactions by the environment dynamical 
correlation effects. 

Let us now examine the effects on the fragment, of the 
average electrostatic field of the rest of the crystal. At 
this point it is important noticing that a majority of the 
strongly correlated materials are ionic or at least partly 
ionic. This is for instance the case of the high- Tc su- 
perconducting cuprates for which the super-conducting 
planes are hole-doped Cu^+Oj" lattices. Similar obser- 
vations can be made for vanadate, manganites, cobaltites 
and even the organic conductors which are charge trans- 
fer salts. It is therefore clear that electrostatic effects 
such as the Madelung potential have to be properly 
treated and that, unlike the field fluctuations, the av- 
erage electrostatic and spin fields of the environment 
strongly affect the fragment under consideration. This 
problem has been carefully studied by Barandiaran and 
Seijo in the case where the fragment borders do not cut 
any covalent bond. They analyzed the crystalline en- 
vironment contributions as Coulomb interactions, quan- 
tum exchange interactions and orthogonality of the frag- 
ment orbitals to the orbitals of the rest of the crys- 
tal (quantum orthogonality interactions) fl5j. Out of 
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these one can separate the short range contributions : 
short range Coulomb repulsions, quantum exchange in- 
teractions and quantum orthogonality, from long range 
ones, that is essentially Coulomb interactions such as 
the Madelung potential. The latter can be easily repro- 
duced using the dominant terms of a multi-polar expan- 
sion, that is for an ionic (or even partly ionic) crystal 
the point charge Madelung potential as described using 
the Ewald [l^ formulation or, as it is more commonly 
done, by a large enough set of point charges modified 
at the set borders by the Evjen's procedure 0. The 
short range contributions necessitate a more sophisti- 
cated treatment since the fragment orbitals should be 
orthogonal to the orbitals (non-included in the calcula- 
tion) supporting the electrons of the rest of the crystal, 
that is essentially with the electrons located on the frag- 
ment first-shell neighboring atoms. The resulting exclu- 
sion of regions in the K'^ space for the fragment electrons 
wave-function has proved to be well treated by total ion 
pseudo-potentials (TIP) approaches. Indeed, the TIPs 
are explicitly derived in order to render on the fragment 
orbitals the quantum orthogonality, quantum exchange 
and short range Coulomb potential effects Illlli. As 
for the Madelung potentials, the short range contribu- 
tions acts essentially at the mono-electronic description 
level (Hartree-Fock). Indeed, both the electrostatic po- 
tential and the TIPs modify only the mono-electronic 
integrals. Thus these effects affect the Fock operator and 
strongly influence the shape of the fragment orbitals and 
their energies. 

While the determination of the TIPs modeling the 
short range contributions is quite unambiguous, the de- 
termination of the Madelung field is more tricky. Even- 
though experimental determinations of the electrostatic 
potential at points locations may be available, this is usu- 
ally not the case and the question of what are the av- 
erage charges supported by the different atoms remains 
open. Most authors use simple chemical considerations 
such as the standard atomic oxidation states and charge 
equilibrium to derive the formal charges supported by 
the different atoms in the crystal. For instance, Nd^^, 
Ou?^ and O^^ in the neodymium-doped superconduc- 
tor Nd20uOi, or Na^+, ^^-5+, 02- ^^ic famous in- 
organic spin-Peierls a'NaV205 compound. Such simple 
considerations assume totally ionic systems. However 
in most systems the ionicity is not as strictly total and 
for systems presenting iono-covalent bonds the field cre- 
ated by the formal charges does not accurately repro- 
duce the real Madelung potential. The problem of deter- 
mining the degree of polarization of such systems or the 
real charges supported by the atoms is therefore crucial. 
That was for instance the case in the a'NaV205 com- 
pound which was commonly assumed in the literature to 
be Att^+yj^'^^Og", while crystal ab initio mean-field cal- 
culations (Hartree-Fock calculations with valence triple 
zeta -|- polarization type of basis sets (lif run using the 



CRYSTAL code [13) yield Miilliken atomic charges of 
^^0.95+^^2.3+ q0.95-q1.3-qi. 15-^ It is well known that 

the Miilliken charges are strongly dependent of the cho- 
sen basis set and in particular of its spatial extension. 
However, the above mean-field calculations clearly ex- 
hibit a strongly-covalent triple bond between each vana- 
dium atom and its corresponding apical oxygen ^^li & 
strong donation from the apical oxygen to the vanadium 
atom — in order to form the third bond — and a weaker 
donation of the other oxygen anions to this V = O sys- 
tem. This result, which couldn't have been derived from 
simple chemical considerations, have been latter con- 
firmed by embedded fragments Complete Active Space 
Self Consistent Field (CASSCF) calculations 0. As 
a conclusion, the choice of the atomic charges generat- 
ing a physically correct Madelung field is a crucial issue 
that cannot be solved by simple charge transfer consid- 
erations. 

The problem can be solved by going back to the ques- 
tion of how the electrostatic potential acts on the frag- 
ment wave- function. As already mentioned, the poten- 
tial generated by the rest of the crystal acts through the 
mono-electronic integrals, and thus, on the orbital shapes 
and orbital energies through the Fock operator definition. 
All these effects belong to zeroth-order description of the 
fragment and are present at the Hartree-Fock level. It 
is therefore possible to calibrate these effects by a di- 
rect comparison of the Hartree-Fock description of the 
embedded fragment and the Hartree-Fock description of 
the infinite crystal. Indeed, the density of states (DOS) 
of the infinite crystal {p'^'^'^''{£)) can be projected on the 
atomic orbitals involved in the fragment under considera- 
tion and compared with the embedded fragment orbitals 
atomic contents and energies. It is easy to define a pro- 
jected density of states (PDOS) function for the finite 
fragment in a similar way as in an infinite system. The 
major difference between the two is that the finite system 
PDOS is a discrete function. If p^J^y^iE) is the PDOS of 
the entire crystal on the atomic orbital /i. 




{k)5{£~E,,{k))dk (1) 



the embedded fragment density of states, pj/'°'^{£), can 
be defined in a similar way as 

= E E A- ^''M ^(^ - (2) 

i I' 

where BZ refers to the Brillouin zone, Vbz to its volume, 
k is the wave vector, c^n{k) is the coefficient of the n**^ - 
band crystalline orbital at point k on the atomic orbital p 
(the star representing the complex conjugate) and £n{k) 
the associated orbital energy. Finally, /^^ is the usual 
coefficient of the j**^ molecular orbital of the fragment on 
the atomic orbital p and £j the associated orbital energy. 
If the embedding correctly reproduces the average effects 
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of the rest of the crystal on the fragment, the discrete 
function pjj''^^{£) should be located at the same energetic 
positions as the crystal PDOS p^''^^(f), and the main 
contributions of both functions should coincide. 



SOME PRACTICAL EXAMPLES 

As already mentioned, the most usual practice used to 
define a fragment embedding consists in 

• building the Madelung potential using formal 
charges, 

• using TIPs on the first shell of cations in order to 
set up the main exclusion effects. 

In the following subsections, we will analyze on different 
examples the validity and limitations of such a procedure, 
as well as the influence of the different parameters on the 
embedding quality. 



The importance of a covalent bond 
compound 



the NaV205 



The a'NaV205 crystal is formed by layers of VO5 
square-pyramids stacked along the c axis. The oxygen 
atoms of the pyramid basis form a quasi-square planar 
lattice along the a and b directions. The pyramids are 
alternatively pointing on top and below these planes. Pe- 
riodical vacancies of the VO top of the pyramids are 
replaced by Na'^ ions forming chains parallel to the b 
axis (see fig. ^1. This system can be modeled as planes 
of weakly coupled ladders, the low energy physics being 
supported by the d^y magnetic orbitals of the vanadium 
atoms and by the py orbitals of the bridging oxygen atom 
on the ladders rungs [l^ . 

a ■ ■ a 

FIG. 1: Schematic structure of a NaV205, (a) along the (a,c) 
plane and (b) along the (a,c) plane. The oxygen atoms are 
denoted by open circles, the vanadium atoms by filled circles 
and the sodium atoms by dots. 

Following a simple chemical analysis, the system was 
given for a long time as a fully ionic system, with a for- 
mal charges repartition of Na^^ (¥2)^^ 0'^~ and a mixed 
valency on the vanadium atom. Later on, quantum 
chemical calculations Tl] showed that the vanadium- 
apical oxygen {V — Oap) bond is in fact a multiple 




strongly covalent bond with a small degree of polar- 
ization (about l/2e). The system should then be seen 
as {V = OapY^ ™ molecules on an oxygen O^^ 
square lattice. A modification of the charge reparti- 
tion in the average iVOapY'^^ group, from y^-^+O^^ 
to V^^O^p^ is considerable and one can expect that it 
strongly affects the potential felt by the magnetic or- 
bitals. Indeed the electrostatic potential (computed by 
the Ewald summation procedure) felt by the atoms sup- 
porting the active orbital goes from PotiV) = 1.735a.u. 
and Pot{Obr) = — 1.103a. u. (where Obr refers to the oxy- 
gen atom bridging the two vanadium atoms on the rungs) 
when formal charges are used to Pot{V) = 1.242a.u. 
and Pot(Obr) = — 0.930a. u. when the polarity of the 
Vanadium-apical oxygen bond has been corrected. One 
notices immediately that the difference of the electro- 
static potential felt by the V and Obr is increased by 
more than 30% when formal charges are used and thus 
that the energy difference between the Vdxy orbital and 
the ObrPy orbital should be strongly affected. 

Figure 13 reports the projected densities of states on 
the rung active orbitals as computed at the Hartree-Fock 
level from a periodic calculation, and from three fragment 
calculations with different embeddings : without embed- 
ding (isolated fragment), with an embedding using the 
formal charges Madelung potential and with an embed- 
ding using the charges corrected to take into account the 
covalency of the V = Oap bond. The rung fragment has 
been chosen as including the VOV rung and the set of 
first neighbors atoms (see fig. . One sees immediately 





FIG. 2: Crystal fragment used for the calculations of the rung 
parameters, along the (a,b) and (a,c) planes. The open circles 
represent the oxygen atoms and the filled circles the vanadium 
atoms. 

that the isolated fragment DOS are very different from 
the periodic DOS, emphasizing, if necessary, the crucial 
importance of the embedding. Indeed, in the isolated 
fragment, even the Fermi level orbitals — Highest Oc- 
cupied Molecular Orbital (HOMO) and Lowest Unoccu- 
pied Molecular Orbital (LUMO) — are not the correct 
ones, i.e. the bonding, anti-bonding combinations of the 
vanadium d^y orbitals. The fragment, embedded with 
the formal charges Madelung potential, yields a reason- 
able representation of the occupied orbitals, the HOMO 
is correctly of vanadium dxy nature and the bridging oxy- 
gen Py orbital is at the correct energy difference. How- 
ever, the unoccupied part of the spectrum is much more 
troublesome, in particular the contribution to the virtual 
orbitals of the vanadium d^y orbitals is totally unreal- 
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istic. We will see in table |I] how this point can affect 
the local spectroscopy of the system and in particular 
the local-singlet local-triplet energy difference associated 
with the super-exchange coupling. Finally the embed- 
ding where the polarity of the V = Oap bond is set to a 
correct value yields reasonable accuracy of the DOS for 
both the occupied and unoccupied parts of the spectrum. 



The influence of geometrical details 
incommensurate Sri4Cu2404i 
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FIG. 3: Projected DOS on the magnetic atomic orbitals of 
the vanadium atoms and the bridging oxygen atom in the 
rung geometry. From top to bottom : periodic crystal HF 
calculation, isolated fragment, embedded fragment with for- 
mal charges and embedded fragment with the correct polarity 
of the V = Oap bond. The dotted vertical line displays the 
Fermi level. 



Table m displays the results for the local-singlet local- 
triplet excitation energy (J) and the first local doublet- 
doublet excitation energy {2t) of the rung fragment, sup- 
porting respectively 4 and 3 magnetic electrons on 3 
magnetic orbitals. One sees immediately the crucial role 
played by the embedding, especially on the effective mag- 
netic exchange integral. Indeed the singlet-triplet excita- 
tion energy is increased by more than 70% when formal 
charges are used. 



In the light of the previous example one may won- 
der how much the details of the Madelung part in the 
embedding can influence the computed local excitations. 
This point is of particular importance for the incommen- 
surate systems composed of several electronically quasi- 
independent subsystems with incommensurate periodic- 
ity. Indeed, it is of crucial importance to know up to 
which point the incommensurability influences the peri- 
odicity of the local effective interactions within one of the 
subsystems. In other words, what is the amplitude of the 
magnetic exchange, hopping, etc. modulation in one sub- 
system due to the electrostatic influence of the other sub- 
system. An interesting candidates to study this problem 
is the super-conducing ladder-chain copper-oxide family 
Sr xCaxCu24:0 41. Indeed, these compounds are consti- 
tuted of alternate layers of weakly-coupled doped spin- 
ladders and weakly-coupled doped spin-chains (see fig- 
ure El- We will use for this study the totally saturated 
X = compound in an ideal geometry where the mod- 
ulation of the chains (resp. ladders) geometry around 
their average structure has been ignored. In this ideal 
system the two isolated subsystems are perfectly period- 
ical and the modulation of the local excitations energies 
along the chain (resp. the ladder) can only be attributed 
to the non-periodicity of the electrostatic infiuence of the 
other subsystem. 
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Charge repartition 
in the Madelung Potential 


J (meV) 2t (meV) 


Na'+ (V*-5+)2(02;)2(0^-)3 


-510 -1162 
-294 -1076 



TABLE L Local excitations energies on the rungs of 
the Qf'A'aVbOs compound as a function of the embedding 
Madelung potential. J singlet-triplet excitation energy of the 
embedded (VO)2^(0^~)7 fragment, 2t first double-doublet 
excitation energy of the embedded {VO)^'^ {0^~)7 fragment. 



FIG. 4: Schematic structure of the incommensurate 
5^^1401(24041 compound along the (a,c) and (b,c) planes. Note 
that the ladder and chain planes have been superposed in the 
(a,c) view. 

The incommensurability of such a system does not al- 
low us to run benchmark periodic calculations in order to 
check the validity of our embedding. However in the cop- 
per oxides, the oxidation numbers of the different atoms 
{Cu, O as well as alkaline-earth elements) do not seem to 
raise much controversy. Indeed, both the chemical analy- 
sis and the ab initio calculations on similar systems such 
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as the chain system Sr20uO^ [i^, the ladder systems 
SrCu20^ li^, or even the perovskites super-conducting 
parent compounds such as the HgBa20un02n+2+s fam- 
ily ^> yield Qi^+ and oxidation states, which we 
will assume for the present system. 

Figure shows the modulation of the local singlet- 
triplet energy difference (computed at the Complete Ac- 
tive Space Self Consistent Field (CASSCF) level) on the 
chain subsystem fragments as a function of their ref- 
erence unit cell along the chain subsystem. The CAS 
space includes either only the copper magnetic d^2_y2 
orbitals (the x and y axes corresponding to a and c crys- 
tallographic axes in figure or as well the orbitals of 
the bridging oxygen atoms — supporting the through- 
bridge, super-exchange mechanism. As for the vanadium 
oxide the embedded fragments have been chosen such 
as to include two nearest neighbors magnetic sites (here 
the copper atoms), the bridging atoms (here one oxygen 
atom on the ladders and two oxygen atoms on the chains) 
and the first neighbor atoms to the previously cited ones. 
The ladder and chain fragments are and Cu20q 

respectively. Unlike the change in the polarization of 
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FIG. 5: Effective magnetic exchange integrals between chain 
nearest neighbor copper atoms, computed at the CASSCF 
level. Circles : the CAS only contains the magnetic cop- 
per d^2_y2 orbitals, stars : the CAS contains in addition 
the two bridging p oxygen orbitals. The shift between the 
chain and ladder unit cells reference point corresponds to 
7n/f0 — Int [n7/10], where n is the present graph abscissa. 

the V = Oap bond in a' NaV205, the modification of the 
Madelung field resulting from the relative shift between 
the ladders and the chains positions induces a very small 
(negligible) effect on the local excitation energies. In- 
deed, the maximal relative variation of the singlet-triplet 
energy difference is smaller than 1.5% along the chain, 
and negligible along the ladders (not shown). This point 
can be easily understood by looking at the Madelung 
potential felt by the different copper atoms along the 
chain or the ladder. Indeed, the potential modulation 
is reported in table ^ One sees immediately, that the 
Madelung potential variation induced by the other sub- 
system is very small, being of the order of 10~*a.u. along 
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TABLE 11: Electrostatic potential in atomic units for the cop- 
per atoms along a chain and one side of a ladder. The in- 
commensurate structure have been approached by the usual 
10 : 7 commensurate one in the chain/ladder direction. The 
10 chain cells, 7 ladder cells mismatch is of 0.145A out of 
27.45A. 

the chains and of 10~^a.u. along the ladders. However, 
it should be noticed that the relative variation of the 
spectroscopic parameters, such as the effective magnetic 
exchange, is two orders of magnitude larger than the rel- 
ative variation of the Madelung potential, pointing out 
the extreme sensitivity of these local excitation energies 
to the electrostatic potential. 

The influence of anions TIPs : the Pt2l[dta)4, 

Finally we will analyze the importance of the short- 
range effects, as modeled by the TIPs, on local excita- 
tions. These effects, excluding the fragment electrons 
from the space supposed to be occupied by the elec- 
trons of the rest of the crystal, are much more impor- 
tant for cations than for anions. Indeed, while the pos- 
itive charges of the former have a tendency to attract 
the fragment electrons, the negative charges of the lat- 
ter already repulse them. Therefore, we will not discuss 
the cations case for which it is commonly admitted that 
these short range exclusion effects are crucial, and we will 
rather analyze the effects due to anions, more subject to 
controversy. 

For this purpose we have chosen the Pt2l{dta)i com- 
pound l23| since the influence of the large /~ ion, first 
neig hbor of the magnetic centers, can be expected 
to be rather important. The Pt2l{dta)i system is a 
quasi-unidimensional system built from the alternation 
of platinum dimers and iodine anions. The coordination 
spheres of the Ft atoms are completed by four /i-bridging 
dithioacetate {dta) ligands (see figure IHJl. The iodine an- 
ion is located at the mid center position between two 
Pt2 dimers. Its distance to the closest platinum atom is 
only 2.98A, while the sum of its ionic radius (2.16A) and 
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the ionic radius of the R'^+ cation (0.98A) is 3.14A. The 




FIG. 6: Pt2l{dta)4 crystal structure. The large empty circles 
are the Pt atoms, the large dark gray circles are the / atoms, 
the middle size light gray atoms are the 5* atoms, the small 
black circles are the C atoms and the tiny black circles are 
the H atoms. 

computed fragment is composed of two Pt atoms bridged 
by the four dta ligands. In this system the electronic 
structure is supposed to be essentially supported by the 
rfz2_|(^2+y2) orbitals of the Pt atoms (the quantification 
axis being set to be the chain axis). The local singlet to 
triplet excitation energy of the [Pt2(dta),^'i^ system has 
been computed with and without TIP [2J| on the first 
iodine neighbors on both sides of the fragment. The sys- 
tem has been found to have a metal-ligand triplet ground 
state. It is noticeable that the singlet, first excited state is 
also a metal-ligand strongly open-shell state. The triplet 
to singlet excitation energy, computed at the CASSCF 
level (the CAS being defined by 4 electrons in 3 orbitals) 
has been found to be of 46meV when the /~ anions are 
only represented by a repulsive potential, that is by 
punctual charges. When TIPs are used, it becomes 
50meV. As can be expected, the TIPs destabilize more 
the singlet state than the triplet state. Indeed, while 
the two states have the same average number of elec- 
trons, the larger charge fluctuations in the singlet state 
increases the probability of the configurations having an 
extra electron on one of the Pt sites, yielding larger repul- 
sive contributions due to the /~ TIPs. It can therefore be 
expected that on systems where the singlet state is less 
open-shell than this one, the first neighbor anions exclu- 
sion effects are larger than in the present case, reducing 
the singlet to triplet excitation energy when the coupling 
is anti-ferromagnetic and increasing it for ferromagnetic 
couplings. However, the absolute value of the TIPs effect 
(4meV) is quite small and can be expected never to be 
very large (since in the present case the iodine-platinum 
distance is already smaller than the sum of the respective 
ionic radii). One can therefore conclude, that the cations 



TIPs will be of importance essentially when the excita- 
tion energies sought at are small (of the order of the tens 
of milli-electron- Volts or smaller). 

CONCLUSION 

We have studied the effects of the environment on the 
local spectroscopy of strongly correlated extended sys- 
tems, as treated in embedded fragment methods. A par- 
ticular interest have been devoted to strongly correlated 
systems and the local effective interactions (magnetic ex- 
change, hopping, repulsion, etc) between the active elec- 
trons responsible for the low energy physical properties. 
We have analyzed the different contributions of the rest of 
the system on the considered fragment, originating both 
from the average electrostatic and spin fields and from 
theirs fluctuations. A multi-polar analysis showed that 
the dominant terms of fields fluctuations contributions 
scales as a charge-dipole contribution (1/i?^). These 
terms decrease rapidly as a function of the fragment re- 
sponse excitation distance R and can therefore be ne- 
glected except from the short range contributions. On 
the contrary, the effects of the average fields are quite 
strong and their long range contributions, originating 
in the Madelung field decrease very slowly. These ef- 
fects (both the short ranged quantum orthogonality and 
quantum exchange and the long range Madelung poten- 
tial) can however be treated at the mono-electronic level 
and thus act essentially on the orbitals shapes and en- 
ergies. We have therefore defined a discrete function, 
equivalent for the finite systems to the Density of States 
in periodic materials, and based a criterion for testing 
the embedding quality on the concordance of the infinite 
system DOS and its discrete equivalent on the embed- 
ded fragment. A detailed analysis of different systems 
have shown that unlike what is commonly assumed in 
the community, the definition of the Madelung potential 
by the formal charges issued from simple chemical con- 
siderations is not accurate. Indeed, these charges usually 
overestimate the charge transfers, often miss the possibil- 
ity of partial covalent bonds, and in any cases are unable 
to evaluate the degree of polarization of the latter. From 
the study of the NaV205 compound, we have shown that 
the crucial importance of these problems. Indeed, the 
orbital shape and energy (and in consequence to the dis- 
crete density of states) are very sensitive to these charge 
transfer problems, as well as the local excitation energies. 
Let us point out that the singlet to triplet local excita- 
tion energy on the NaV20z rungs differ from a factor 
1.7 when the Madelung potential is defined from formal 
charges and when the partial covalency of the vanadium 
apical-oxygen bond is correctly defined. On the contrary, 
other effects such as the short range exclusion effects due 
to the fragment first neighbor anions are very small and 
dominated again by the \/ R electrostatic repulsion term. 
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In conclusion wc would like to emphasize again the 
importance of a correct evaluation of the charge transfer 
processes and a validation of the embedding by check- 
ing the embedded fragment orbital shapes and energies 
against the infinite system density of states. 

We thank Dr. D. Maynau for providing us with the 
CASDI set of programs used in some of the presented 
calculations. 
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